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Chiral Symmetry Breaking in the 3-d Thirring model for small Nf. 
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a Centre de Physique Theorique, CNRS Luminy, F-13288 Marseille Ccdex 9, France 

We study the dynamical breaking of chiral symmetry in the 3-d Thirring model for a small number of fermion 
species. The critical point is identified by fitting lattice data to an equation of state. The spectrum of the theory 
is studied to confirm the phase structure of the model. 
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1. Introduction 

The Thirring model is a three dimensional QFT 
with Nf fermion species and a current-current 
four fermion interaction. As such, it is obviously 
non-renormalisable in the usual perturbative ex- 
pansion. It has been pointed out several years 
ago [0 that a continuum limit can be defined 
for some theories in a non-perturbative (NP) ap- 
proach. For the Gross-Neveu model, a four-fermi 
theory with a (pseudo-)scalar interaction, the ex- 
istence of such an UV fixed point has been shown 
within the 1 /Nf expansion and confirmed by nu- 
merical simulations. The critical point separates 
a strong coupling phase with chiral symmetry 
breaking from a symmetric weak-coupling one. 

For the Thirring model, solid conclusions can 
not be drawn from analytic results as different 
NP techniques yield different results. We aim to 
obtain a complete description of the phase struc- 
ture of this model by lattice simulation, using 
staggered fermions and a HMC algorithm. The 
motivations, the lattice formulation of the the- 
ory and the technical details of our analysis were 
presented in great details in earlier works to 
which we refer the interested reader. This con- 
tribution summarizes our most recent results for 
Nf = 4,6. 

The order parameter for detecting chiral sym- 
metry breaking is the chiral condensate (xx) • Its 
behaviour around the critical point is described 
by fitting to an equation of state (EOS) jj). 
Finite-size effects are included in the EOS and 
the critical exponents are obtained for Nf = 4. 



The results are compared with those for Nf = 2. 
Once again, lattice data show a clear signal of 
chiral symmetry breaking and fit rather well the 
RG predictions. 

Moreover the chiral symmetry breaking transi- 
tion should manifest itself in the behaviour of the 
susceptibilities and in the spectrum of the theory. 
This issue was already addressed for Nf = 2 in Q . 
Preliminary results for Nf — 4 are reported here. 
The general pattern, which had already emerged 
in our previous works, is confirmed by the new 
results. We briefly outline the next developments 
in our program. 

2. Fitting to the Equation of State 

The EOS is a functional relation between the 
order parameter, the coupling, <?, and the bare 
mass, m, in the regularised theory, which can be 
derived form the RG equations in the vicinity of 
a fixed point (see e.g. H)). It therefore provides a 
relation to which lattice data can be fitted. The 
results of this approach for Nf = 2 and fixed 
lattice size were presented in || . In order to take 
into account finite size effects, it is sufficient to 
consider the inverse linear size of the lattice, l/L 
as an additional scaling field of dimension 1 . The 
exact functional form of the EOS is unknown, but 
it can be expanded in Taylor series, yielding a 
fitting function for the lattice data with only a 
few unknown parameters: 



m = B(xx) 5 + A(t + CL~^)(XX) 5 - 1/0 



(1) 
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where t = 1/g 2 — \/g\. 

For Nf — 4, we fit to the EOS data from lattice 
whose sizes range from 8 4 to I6 4 . The results are 
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Table 1 



Results from fits including finite 


size scaling. 








Parameter 


Fit 




Parameter 


Fit 


Nf = 2 1/gi 


1.92(2) 


iV/ = 4 


Wc 


0.70(2) 


S 


2.75(9) 




s 


3.41(31) 


P 


0.57(2) 




P 


0.41(10) 


n 


0.60(2) 




n 


0.34(20) 


V 


0.71(4) 




V 


0.61(10) 


A 


0.334(7) 




A 


0.70(5) 


B 


2.7(3) 




B 


5.8(2.1) 


C 


2.1(7) 




C 


2.9(2.2) 


X 2 /d.o.f 


1.76 




X 2 /d.o.f 


0.80 



summarized in Tab. [l]. 

It is interesting to stress that the results from 
the fit do take into account the effects due to the 
finite volume of the lattice and should therefore 
be interpreted as values in the thermodymamical 
limit. The values for the critical coupling and the 
critical exponents are stable with respect to the 
previous results obtained for Nf = 4 at fixed lat- 
tice size 0. Such a stability in our fits can be 
seen as a sign that the lattice data really follow 
the behaviour expected in the vicinity of an UV 
fixed point. The critical exponents differ signifi- 
cantly between Nf = 2 and Nf — 4, confirming 
a result that had already emerged from the fixed 
size fits. We should be able to reduce the er- 
rors on the fitted values as statistics increase over 
the next months. It is worthwhile to stress once 
again that, in order to identify unambiguously a 
fixed point, we should study the renormalised tra- 
jectories in coupling space, for which we have to 
measure 3- and 4-point correlation functions. 

3. Susceptibilities and spectrum 

Independent information on the phase struc- 
ture of the theory can be obtained from the sus- 
ceptibilities and the spectrum of the theory. Once 
again precise definitions and the results for the 
Nf — 2 case are given in previous works and 
we want to concentrate on the new results for 
Nf = 4, 6. We briefly recall the definitions and 
the expected behaviour for some quantities of in- 
terest. 

The transverse and longitudinal susceptibilities 



are defined respectively as: 

XI = \ E<XX(0) Xx(aO> ( 2 ) 

X 

Xt = ^V^XexlO) X£x{x)) c = — (xx) (3) 
V '-^ m 

X 

They are related to the scalar and pion masses: 
Xi,t = Z^/M* (4) 



Another interesting quantity is their ratio R = 
Xi/xt- If there is a critical point in the phase dia- 
gram where chiral symmetry is dynamically bro- 
ken, then: 



lim R 

m—>0 




g 2 > g 2 c- 
g 2 < g'c- 



(5) 



The results for the susceptibilities are shown in 
Fig. 0. Even though the results were produced 
at non- vanishing bare mass, i.e. in presence of 
an explicit symmetry breaking term, the pattern 
reflects the existence of a phase transition. Once 
again the near degeneracy of \i an d Xt i n the 
symmetric phase (t > 0) supports the hypothesis 
d — = 1, used in the fits reported in Tab. [[]. 

In the strong coupling regime (t < 0), the 
transverse susceptibility is larger than the lon- 
gitudinal, reflecting the fact that the pion is al- 
most a Goldstone boson in that regime. The same 
trend is confirmed by the masses extracted from 
the aymptotic behaviour of two-point functions. 
Preliminary results are reported in Fig. || the er- 
rors on the scalar mass is still very large and has 
not been included in the plot. A more detailed 
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Figure 1. Transverse (circles) and longitudinal 
(diamonds) susceptibilities for Nf = 4, to — 0.01, 
on a 16 3 lattice. The data are plotted against a 
rescaled reduced coupling t = g 2 (l/g 2 — l/fic)- 



analysis, based on larger statistical samples, is 
currently in progress and we hope to report soon 
on the details of the spectrum for Nf — 4. It 
is nonetheless worth to remark that the general 
pattern is consistent with the existence of an UV 
fixed point. In particular, we see no sign for the 
unconventional spectrum in the symmetric phase 
associated to conformal phase transitions, as sug- 
gested in [El. 
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Figure 2. Fermion (circles), pion (diamonds) and 
scalar (triangles) masses vs. 1/g 2 , for Nf = 4, 
to = 0.01, on a 16 3 lattice. 



In order to distinguish between different sce- 
narios of chiral symmetry breaking (correspond- 
ing, e.g., to different solutions of the Schwinger- 



Dyson equations), several values of Nf need to 
be studied. New preliminary results for the chi- 
ral condensate are reported in Fig. [| The lattice 
used is larger than in our earlier works Q and the 
behaviour of the condensate as the bare mass to 
is decreased suggests the existence of a symmetry 
broken phase at strong coupling. We will report 
more systematic results in the near future. 



Figure 3. Chiral condensate vs. 1/g 2 for Nf = 6, 
m = 0.01 (circles) on a 16 3 lattice and m — 0.02 
(squares) on a 12 3 lattice. 
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